Control of Local Relaxation Behavior in Closed Bipartite Quantum Systems 
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Institut fur Theoretische Physik 1, Universitat Stuttgart 

We investigate the decoherence of a spin 1/2 subsystem weakly coupled to an environment of many 
spins 1/2 with and without mutual coupling. The total system is closed, its state is pure and evolves 
under Schrodinger dynamics. Nevertheless, the considered spin typically reaches a quasi-stationary 
equilibrium state. 

Here we show that this state depends strongly on the coupling to the environment on the one hand 
and on the coupling within the environmental spins on the other. In particular we focus on spin 
star and spin ring-star geometries to investigate the effect of intra-environmental coupling on the 
central spin. By changing the spectrum of the environment its effect as a bath on the central spin 
is changed also and may even be adjustable to some degree. We find that the relaxation behavior 
is related to the distribution of the energy eigenstates of the total system. 

For each of these relaxation modes there is a dual mode for which the resulting subsystem ap- 
proaches an inverted state occupation probability (negative temperature). 

PACS numbers: 05.30.-d, 03.65.Yz, 75.10.Jm 



I. INTRODUCTION 

In a composite but closed quantum system in which 
a smaller "central system" is weakly coupled to a larger 
"environment" , most of the (pure) states of the total sys- 
tem for a given energy (and possibly some additional con- 
straints) exhibit properties of thermal equilibrium states 
with respect to the smaller part 0,0 j he. there exists a 
so-called dominant region in Hilbert space in which the 
entropy of the central system is close to its maximum 
value under the given constraints. For a completely "un- 
structured" coupling to the environment, the state of the 
central system, starting away from equilibrium, shows 
decoherence and thermalization and, ultimately, a quasi- 
stationary equilibrium situation is reached which is de- 
termined only by the spectrum of the environment [j| . 

In this paper we investigate to what extent structured 
systems will also exhibit this thermalization behavior. 
We focus, in particular, on a single central spin 1/2 par- 
ticle coupled to a relatively small environment of spin 1/2 
particles. Recently, the properties of spin systems of dif- 
ferent structure (rings, stars and others) have been sub- 
ject of extensive interest. A lot of work has been done 
on the question of entanglement 0, 0, pi p i I S |jj , their 
relaxation behavior has been addressed \WL [jjj and vari- 
ous techniques were suggest ed to make any spin interact 
with any other spin [Ho , Il3| . 

Here we show that by choosing different types of cou- 
pling between the central system and the environment on 
the one hand and within the environment on the other 
the equilibrium state finally reached can be controlled. 
There is, in particular, a qualitative difference between 
the relaxation dynamics and the equilibrium of the cen- 
tral spin for interacting and noninteracting environments. 
We relate this equilibrium to the spectral structure of the 
total system. 



II. SPECTRAL TEMPERATURE 

Consider a weakly coupled bipartite quantum system 
consisting of a small system S and a large "container" or 
environment C with the Hamiltonian 



H = H 



aH" 



(aH int ) « (H s ) , (H c ) . (1) 



The behavior of the system is completely described by a 
Hilbert space vector \i/j(t)) and its evolution under the 
Schrodinger equation. Under weak constraints, the state 
of S after relaxation is determined by Jaynes' principle, 
taking energy as the only relevant observable 0, 0] • The 
temperature of the respective canonical state can be pre- 
dicted from the degeneracy structure of the environment 
alone. 

Fig. ^a) shows a two-level system in resonant con- 
tact with an environment consisting of two "energy 
bands" fc, k' of degeneracies gu and g^, respectively, 
(gk,gk' S> 1) in a non equilibrium (initial) state. In equi- 
librium, the state of the total system is expected to be 
distributed homogeneously over the whole Hilbert space 
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FIG. 1: Two-level system S in contact with environment C 
consisting of two highly degenerate levels k, k' with degenera- 
cies gk and g^, respectively, (a) Resonance, 6 s — S c , showing 
the state |1) ® \k,m) for which S evolves into a thermal 
state, (b) Slightly detuned energy splitting between S and C: 
Population relaxation in S is prohibited by energy conserva- 
tion, C no longer acts as an effective heat bath. 
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on average, the time-averaged reduced state operator of S 
is thus given by 



1 



P 



(gv |0}(0|+ 3fe |i)(i| 



(2) 



9k + 9k' 

which can be interpreted as a canonical state operator 
with inverse temperature 
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(3) 



In general, the inverse equilibrium temperature imparted 
on S for participating energy levels E of the environ- 
ment C can be defined as the derivative of the logarithm 
of the degeneracy g (E) with respect to energy, 

d\n(g c (E)) 



dE 



(4) 



Note that this is a spectral property only and does not 
imply that the environment was in a thermal state; in 
fact, under the condition of Schrddinger dynamics for 
the total system it is far away from such a state. If the 
degeneracy of environmental states increases exponen- 
tially with energy, s is independent of the initial state 
of the total system (S + C). In this case also systems S of 
higher dimensional Hilbert space will evolve into canon- 
ical states. 

Equations are valid only if energy exchange be- 

tween S and C is possible, i.e. if transitions in both sub- 
systems are in resonance. Under more general conditions 
it is far from clear what numbers should count as the 
pertinent degeneracies gk- One may expect that the ac- 
tion of the environment as an effective heat bath could 
even be turned on and off by tuning the level splitting in 
either system as shown in Fig. ^b). In case of entirely 
detuned transitions the environment would still induce 
decoherence, i.e. phase relaxation in S, thus effectively 
acting as a microcanonical environment. 

At first glance one would expect both parts of the sys- 
tem to be in resonance if the detuning \S S — S c \ was 
smaller than the width of the involved bands in C, i.e. it 
best should be zero. However, if the free Hamiltonian 
of S is renormalized by the environment it is possible 
that the transitions within the central system and the 
environment become somewhat off resonant for vanish- 
ing detuning. In this paper we focus on values of the 
detuning with maximal overlap, i.e. we choose S s — S G 
such that the equilibrium state was as close as possible 
to the expected one from equation J2J). However, for part 
of the following discussion the precise value of the detun- 
ing is not very important, and <5 S could have been chosen 
equal to 5 C . 



III. CHOICE OF ENVIRONMENTAL 
SPECTRUM 

To use the quantum environment as an effective heat 
bath for the central system S with adjustable tempera- 
ture, though, the degeneracy of the environment must 
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FIG. 2: Solid line: Inverse equilibrium temperature /3 s of S 
in units of S c depending on the involved energy band k of 
C. Dashed line: Corresponding population inversion I s in 
equilibrium. The results shown are for N = 50. 



not scale exponentially with energy. Then the equilib- 
rium temperature of the subsystem depends on the ini- 
tial state of the environment via equation J3J and can be 
adjusted by adequately preparing the initial environmen- 
tal state. For the equilibrium state of the system S to be 
sufficiently independent of its initial state, the spectral 
temperature should vary only little over adjacent energy 
levels, though. 

A possible implementation of these requirements is a 
degeneracy varying binomially with energy [l4j |. 
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where we assume resonance of S with adjacent energy 
bands k and k + 1, < k, k + 1 < N. The band 
splitting S c will be taken as a convenient energy scale. 
For N = 50 Fig. [2] shows the equilibrium temperature 
of S as a function of the "band index" k of the envi- 
ronmental initial state and the corresponding population 
inversion I s = (1| p s |1) — (0| p s |0). The inversion of the 
state (01 depends on k and is given by 



I S = (&,) = 



9k - 

gk + gk+i 



(5) 



For sufficiently large N and < k < N, the tem- 
perature T oc 1/ s varies over a wide range of posi- 
tive, infinite and even negative values, although locally, 
i.e. around a given energy band, the deviations from an 
exponential degeneracy-structure are still small, except 
for very small or large k. 



IV. SPIN ENVIRONMENTS 

An environment consisting of N spin 1/2 particles with 
dominant Zeeman-tcrms 
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J Z + H CG , j z = Y, a 
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(mutual interaction H cc within C small) has the desired 
binomial degeneracy, ai^ denotes the Pauli-operator of 
the V th particle. To this environment we couple the cen- 
tral spin with the local Hamiltonian H s = ^-a z . 

We will now consider several different types of inter- 
action H lnt between S and C. To test whether the sys- 
tem shows thermalization for a given H lnt , we study the 
Schrodinger time evolution of the initial product state 
depicted in Fig.^a), 

\*P s+c {t = 0)) = |l) s ® \k,m) c , (6) 

where k denotes an energy band in the environment 
and m indicates one state within this band. If the cou- 
pling is weak and <5 S w S c energy conservation restricts 
the dynamics mainly to a subspace of Hilbert space di- 
mension 



dn = 9k + 9k+i = 




spanned by the states {|l) S ®|fc,m) c , |0) S ®|fc + 1,to') c } 
(to = 1, 2, . . . , <j>fe and to' = 1,2,..., gu+i)- After a suit- 
able thermalization time the reduced state operator of 
the central system p s should be close to the one given 
by ©. 

The subspace of Hilbert space spanned by the gu+gk+i 
state vectors {|1) S ® \k, m) C , |0) S (£> |fc + 1, to') } will be 
called the "accessible subspace" throughout this paper. 
Nevertheless, other states have to be included, if one is in- 
terested in the off-diagonal elements of the reduced state 
of S. The accessible subspace is to be distinguished from 
the subspace that is actually "involved" in the dynamics 
from a given initial state via Schrodinger time evolution. 
This "involved" subspace usually is a subset of the "ac- 
cessible" one. 
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FIG. 3: ^-component of the Bloch vector of S for initial 
state © (N = 14, k — 2, g s = ( 1 3 4 ), g 2 = (")) and ran- 
dom perturbations as a function of time t. The horizontal 
line indicates the expected mean value I s = (ct z ) = -3/5 
(see Q). 

455. Regardless of the small number of dimensions one 
can observe relaxation to the expected mean value indi- 
cated by the horizontal line. The average over all times 

W) = lim i f ' Tr{<7 2/ 3 s (i)}dt 

T^oo 1 J _ T / 2 



is (a z ) = —0.598 for the simulation of Fig-El Considering 
the small number of spins this is remarkably, but not un- 
expectedly, close to the expected average —3/5 = —0.6. 

A change of the coupling strength a alone affects only 
the timescale of the dynamics, provided the coupling 
stays weak, since all energy eigenstates of the uncoupled 
system in the accessible subspace are degenerate. The 
detuning can be changed slightly without disturbing the 
system qualitatively as long as the introduced level split- 
ting is considerably smaller than the interaction energy. 



A. Random Interaction 

To test the expectations of the previous sections 
we first take H cc = and H mt as a hermitian 
random matrix form the Gaussian unitary ensemble 
in 2 JV+1 -dimensional Hilbert space with the distribu- 
tion w(H int ) = Cexp(-ATr{H int } 2 ) [l5|. It has been 
shown previously Q that an initial product state evolves 
into the desired canonical state under the time evolution 
generated by the given Hamiltonian. This is not surpris- 
ing, since, in general, the complete information about the 
environment being constructed from spin 1/2 particles 
is lost and energy remains the only conserved quantity. 
Most states in the accessible region of Hilbert space are 
equilibrium states with respect to the central system. 

Fig. |3 shows the time evolution of the z-compo- 
nent (<x z ) = Tr{a z p s (t)} of the Bloch vector for N = 
14 particles in C and initial state © with k = 2, S s = S c 
and very weak coupling a = <5 c /5000. The accessible 
subspace of Hilbert space is of dimension c?h = gi + 53 = 



B. Spin-Star Configuration 

The situation is quite different if the interaction is 
structured. Neglecting intra-environmental interaction 
(i.e. H cc = 0), the most general Hamiltonian for cou- 
pling the central spin to each spin in the environment is 

3 N 

Even highly symmetric realizations of this class like the 
Heisenberg XX interaction ^(| show dissipation and 
decoherence with respect to S for highly mixed initial 
states [lOj. 

However, this is no longer true for the pure initial 
state ©. Fig. 0] shows the time evolution of (a z ) of S 
for N = 14 particles in the environment and 6 « 6 C = 
5000a. The coefficients 7^. have been chosen randomly 
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FIG. 4: 2-component of the Bloch vector of S for initial 
state © (N = 14, = 2) for two different choices of # int : 
Random two-body interactions, 5 s = 1. 000825° (solid line), 
S s = 0.9992(5° (dashed line). 



FIG. 5: z-component of the Bloch vector for initial state © 
(N = 14, k = 2) with stronger (solid line), weaker (dashed- 
dotted line) and without (dashed line) intra-environmental 
coupling. See the text for the respective values of 7. 



from a normalized Gaussian distribution. Here oscilla- 
tions with a significantly larger amplitude than in the 
previous section are observed because a smaller fraction 
of Hilbert space is involved in the time evolution. Also 
the time- averaged (a z ) differs considerably from the ex- 
pected value indicated by the horizontal line. This av- 
erage over all times is (<r z ) = —0.0976 for the dashed 
and (a z ) = —0.261 for the solid line. 

As stated earlier, 5 s — S c was chosen to "optimize" 
the relaxation behavior with respect to approximating 
state (J2J. Any other value than that used would lead to 
an average reduced state of S even further away from the 
expected one. 

C. Ring-Star Configuration 

Now we extend the spin-star model by introducing mu- 
tual interaction of the environmental spins. Next neigh- 
bor coupling in the form of the quantum fsing chain, 

N 

is added to the system discussed in the last section 
(assuming periodic boundary conditions, i.e. = 
a£ '). Fig. [3] compares the time evolution of the spin- 
star model as discussed in the last section (dashed line, 
same as dashed line in Fig. [3J with the same system to 
which this H cc has been added (solid line with 7 = 3ev 
and dashed-dotted line with 7 = 3a/2). The oscilla- 
tion amplitude is reduced considerably and the time av- 
era ged v alue is closer to the expec ted o ne. Here we 
get (a z ) = —0.591 for stronger and (a z ) — -0.432 for 
weaker additional intra-environmental interaction. 

The relaxation dynamics and equilibrium state of S 
strongly depend on the coupling within the environment. 
By increasing 7 from zero the amplitude of the oscilla- 
tions is reduced and the average value of (a z ) approaches 
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FIG. 6: (a z ) for different values of 7. Initial state © (N = 14, 
k = 2),S s = 5 c . 

its expected value. This is explicitly shown in Fig. [ft] 
which also shows that the precise value of <5 S — 5 C is not 
too important here. We have found similar behavior for 
XY- and Heisenberg-type interaction. Ising-type inter- 
action (a z ® a z ) is somewhat different: By carefully ad- 
justing the parameters the system approaches an average 
state close to the expected one, yet the fluctuations are 
not decreased considerably. It is important to note that 
interactions involving a z change the effective band split- 
ting in the environment so that the value of the detuning 
becomes significant and has to be adjusted accordingly. 

One may imagine to switch on or off the interac- 
tion H cc or to tune its strength dynamically (via "re- 
focusing" pulses, see e.g. ^3)- ^ n this wa Y the environ- 
ment would become adjustable with respect to its effect 
on the given central spin. This is a phenomenon which in 
"classical" thermostatistics could hardly have been antic- 
ipated. Furthermore, this shows how to control S entirely 
by a specific modification of its environment. 

Finally we note that in addition to showing dissipation 
(relaxation of the diagonal of p s ) superpositions in S are 
decohered (the off-diagonal elements are zero on aver- 
age) , see Fig. [7| 
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FIG. 7: x- and z-components of the Bloch vector for an initial 
superposition (j0) + |l})/\/2 of S as functions of time. Solid 
line: {a x ), dashed line: (<r z ). 



V. SPECTRAL PROPERTIES 

Finally, we try to attribute the various types of relax- 
ation behavior of these bipartite quantum systems to the 
distribution of the energy eigenstates in Hilbert space or, 
more precisely, to the distribution of certain properties of 
the energy eigenstates. Let \e n ) be an energy eigenstate 
of the total system defined by . Then 

e s n = Tr c {\e n ){s n | } 

is the corresponding reduced state of S. Since energy con- 
servation restricts the dynamics to some Hilbert subspace 
as discussed in section llVl (the accessible subspace), we 
only take the energy eigenstates of this subspace into ac- 
count in the following discussion. 

An interesting property of the e„ is the z-component 
of their Bloch vectors X z . n = Tr{cr z e^}. The average 
of all X z .n of a given Hamiltonian in the discussed en- 
ergy band is approximately given by the expected value 
for the equilibrium JSJ . Figures |H1 to ^| show the dis- 
tributions of the X Zjn for each of the models (N = 14 
environmental spins) discussed in the previous section 
over the ensemble of the respective random Hamiltoni- 
ans. These distributions have been calculated numer- 
ically for 100 different randomly chosen realizations of 
each model, i.e. 45500 values in total for 14 environmen- 
tal spins (as in section llV|) . 

Fig. [H] shows the distribution of the values of X z 
for the complete random interaction discussed in sec- 
tion II V Al Because the eigenvectors of these Hamiltoni- 
ans are homogeneously distributed over the unit sphere 
in the accessible subspace, the X z are narrowly peaked 
around —0.6 01- Due to the homogeneous distribution 
of the eigenvectors one expects most of them to con- 
tribute to the dynamics with approximately equal weight 
and therefore (a z ) is expected to be close to the aver- 
age X z . 

For the spin-star configuration shown in Fig.[5]the sit- 
uation is quite different. Although the average over all X z 
is still —0.6, there is a strong peak at X z — —1. All the 
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FIG. 8: Distribution of the z-component of the reduced state 
operator corresponding to the energy-eigenstates for a com- 
pletely random perturbation as discussed in section IIV Al 
14 particles in C. 
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FIG. 9: Distribution of the z-component of the reduced state 
operator corresponding to the energy-eigenstates for the spin- 
star configuration discussed in section llV Bl 14 particles in C. 



energy eigenstates corresponding to the region X z rs — 1 
are approximately product states of the form 

K)«|o) s ®lx) c 

with some state \\) C of the environment (a correspond- 
ing peak for eigenstates of the form \e' n ) « |1) S ® |x) is 
missing here since this state is not part of the accessible 
subspace). For the initial condition given by © these 
states do not contribute to the dynamics, therefore the 
subspace "involved" in the dynamics is only a small frac- 
tion of the "accessible" subspace, leading to larger fluc- 
tuations in the time evolution. Also the average X z of 
the contributing states is larger than —0.6, therefore the 
equilibrium temperature differs from its value for a com- 
pletely random perturbation. The fluctuations could be 
decreased simply by using a larger environment. How- 
ever, the average X z of the contributing energy eigen- 
states, and therefore also the average equilibrium state 
would not be affected. 

The shape of the distribution can be changed by in- 
troducing intra-environmental coupling. For the types 
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FIG. 10: Distribution of the ^-component of the reduced state 
operator corresponding to the energy-eigenstates for the ring- 
star configurations discussed in section HV CI 14 particles in 
C. 

of interaction discussed in the previous section the peak 
is shifted to larger values of X z for increasing coupling 
strength, Fig.llOlshows the result for 7 = 3a. The distri- 
bution is considerably broader than in Fig. |SJ while the 
peak is again close to X z — —0.6 and the "involved" sub- 
space is enlarged. Since no large fraction of the energy 
eigenstates is excluded here by the choice of th e init ial 
state, most of the X z contribute and one expects (a z ) to 
be again close to the average X z . 



VI. "THERMAL DUALITY": RELAXATION TO 
NEGATIVE TEMPERATURE STATES 

For k > N/2 the degeneracy of the environment de- 
creases with energy, dln(g c (E))/dE < 0, and there- 
fore C can act as an environment imparting a negative 
absolute temperature to the central spin. 

We have shown that the environment leads for the con- 
sidered subsystem to Schrodinger relaxation approach- 
ing a quasi-stationary state with positive temperature 
(for k < N/2). These scenarios may be considered as 
generic, if simplified, models for what we typically ob- 
serve in real life. The environments with a finite num- 
ber of states allow also for a different class of relaxation 
modes: As a typical consequence of a quantum approach 
to thermodynamics we find relaxation towards negative 
temperature states as we replace the accessible environ- 
ment band index k by N — k (the relaxation behavior of 
the system is fully symmetrical with respect to the trans- 
formation \s) (g> \k, m) — > \s) ® \N — fc, to), where s means 
logical negation)! Contrary to inversion due to optical 
pumping, these dual states would be dynamically stable 
just as their positive temperature counterparts. Note, 
again, that the environment is not a heat bath of neg- 
ative temperature; in fact, it is far from any canonical 
state. 

Since the population of the energy levels of S at T < 
are inverted (I s > 0), negative temperatures are actu- 
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FIG. 11: z-component of the Bloch vector for the ring-star 
geometry of section II V CI with initial state © (N = 14, k = 
11). 



ally "hotter" than positive ones. Spin systems and ther- 
modynamics at T < have been considered before, see 
e -g- |3E3- I n this regime several statements of thermo- 
dynamics, like the second law, have to be reformulated. 

Here we show relaxation of S to a state of negative 
absolute temperature. The situation is the same as dis- 
cussed in sections IlIII and ITVl but now we use |0 , k, m ) 
and \l s ;k',m c ) with k = 12 = N - 2, k' = 11 = N - 3 
(instead of k — 2, k' = 3) as the accessible subspace for 
the Schrodinger time evolution (see Fig.|2J). 

Fig.llllshows (a z ) (t) for a ring-star geometry of the en- 
vironment as discussed in section llV Cl and demonstrated 
in Fig. [5J As expected, the time averaged inversion is 
close now to 0.6 instead of —0.6. Thus for environments 
consisting of many spins, the regime T < can natu- 
rally be reached for sufficiently high energy in the en- 
vironment. In this way it is possible to study quantum 
thermodynamical effects at negative absolute tempera- 
ture and possibly even heat conduction between embed- 
ded subsystems at temperatures of different sign. This 
"dual" world, though artificial, may contribute to a bet- 
ter understanding of thermal physics based on quantum 
mechanics. 



VII. CONCLUSION 

For the example of a simple subsystem (spin 1/2) 
weakly coupled to a spin environment we have shown that 
the spectral structure of the environment has a strong 
influence on the decoherence of the central system. We 
have compared the decoherence of the central spin due 
to random and structured coupling. 

A noninteracting spin environment does not induce re- 
laxation into the local thermal state expected from the 
band structure, even though energy is the only conserved 
quantity (the spectrum of the total system is non degen- 
erate) . In this case a majority of energy eigenstates do 
not constitute linear combinations of all states permitted 
by energy conservation and therefore the accessible part 
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of Hilbert space is considerably reduced for certain non 
equilibrium initial states like the pure state ©■ 

By introducing additional mutual coupling within the 
environment, thus changing its spectrum, some of the 
properties of the completely random perturbation with 
respect to relaxation of the central system are restored, 
provided the coupling strength is adjusted appropriately. 
By dynamically changing this coupling strength the thcr- 



malizing effect of the environment on the central system 
should be adjustable. 

We have concluded our investigation with some explo- 
rations on thermal duality: Negative temperatures of the 
embedded system appear here as a natural consequence 
of the spectral properties of the environment. 

We thank the Landesstiftung Baden- Wiirttemberg for 
financial support. 
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